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Abstract 

The quantum relativistic Buneman instability is investigated theoretically using a collective 
Klein-Gordon model for the electrons and a cold fluid model for the ions. The growth rate and 
unstable wave spectrum is investigated in different parameter regimes corresponding to various 
degrees of relativistic and quantum effects. The results may be important for streaming instabilities 
involving ion dynamics in very dense plasmas. 

PACS numbers: 52.27.Ny, 52.35.Qz, 67.10.Hk 
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I. INTRODUCTION 



The problem of the stability of electron beams propagating through a plasma is important 
in the context of laboratory beam-plasma experiments jl], inertial fusion schemes (2 ] , in the 
solar corona {sj], astrophysical objects ^,[5], etc. It has been suggested that pulsar glitches 
are due to a streaming instability , where super-fluid neutrons and superconducting 
protons co-exist with relativistic electrons 8]. It was early recognized that electron beams 
propagating through the plasma can give rise to high-frequency electron plasma waves, or low 
frequency ion acoustic waves. The ion motion becomes important if the electrons drift as a 
whole through the plasma and typically becomes unstable if the relative drift speed between 
the electrons and ions is larger than the ion acoustic speed, as in the original Buneman 
instability in an unmagnetized plasma (or electrons streaming along the magnetic field 
lines) and the Farley-Buneman instability in magnetized plasmas [lo| . 

The physics of the linear Buneman instability is fairly well understood in the non- 
quantum (h = 0) and non-relativistic (1/c = 0) cases (h is the reduced Planck's constant 
and c is the speed of light in vacuum). In other situations, relatively thin beams of electrons 
penetrate a plasma consisting of both electrons and ions, giving rise to electron two-stream 



and ion Buneman instabilities in different parameter regimes in the non-relativistic 12| and 



relativistic 



13Nl5| regimes. An electron beam propagating through a plasma can give rise to 



electrostatic two-stream instabilities when the wave vector and electric field are aligned, and 

electromagnetic instabilities obliquely to the beam direction. A typical condition for the ion 

motion to become important in a current-neutral plasma is 7 > aM/(mZi), with 7 being 

the electron beam gamma factor, a the beam to plasma density ratio, Zj the ion charge 

state, and M, m the ion and electron masses [2]. To incorporate quantum effects, collective 

multi-stream Schrodinger 16j| and Klein-Gordon [17| models have been used to investigate 

the electron quantum two-stream instability for non-relativistic and relativistic cases. In the 

non-relativistic regime, the quantum two-stream instability can be physically understood in 

terms of a free energy available due to negative energy modes 181 ] . In addition, low frequency 

linear and nonlinear waves in collisional, non-relativistic quantum plasmas can be derived 

□ 

using multi-stream, carrier-envelope methods 19]. 

In this work, we consider the quantum relativistic Buneman instability, in which the bulk 
electrons stream against the ions, initially at rest. We compare the classical relativistic 
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case with the non-relativistic and quantum well as with the combined quantum 



20| for the electrons 



relativistic case. As a model, we use a collective Klein-Gordon model 
and a cold plasma model for the ions. 

The paper is organized as follows. In Section II, the pertinent dielectric function and 
dispersion relation are deduced, starting from the basic hydrodynamic equations for the 
two-species plasma. In Section III, assuming low frequency and a beam-plasma resonance 
condition, the corresponding Buneman instability is derived, in the case of negligible quan- 
tum effects. Moreover, the instability is intuitively understood in terms of positive and nega- 
tive energy modes. Section IV and V extend the treatment to the quantum non-relativistic 
and quantum relativistic cases, respectively. Section VI contains the final remarks and a 
discussion on the validity conditions of the model. 

II. DIELECTRIC FUNCTION 

The mathematical model is based on a collective Klein-Gordon equation for the electrons 
that is cast in to a set of fluid-like equations, and a cold fluid model for the ions. The 
electrons are described by a Klein-Gordon field ip = Rexp[iS/h], so that R = R(x,t) and 
S = S(x, t) can be viewed as, respectively, the amplitude and phase of the collective electron 



wave function, 
obtained as [l?], 



? rom the Klein-Gordon equation, the evolution equations for R and S are 
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1 c 2 dt 2 dx 2 ) c 2 dt c 2 dt \dt ) dx dx 

where — e and m are the electron charge and mass, c is the speed of light in vacuum, and 
h is Planck's constant divided by 2tt. Equation (JTJ) plays the role of a quantum relativistic 
electron continuity equation, while Eq. (J2J) is a relativistic Hamilton- Jacobi equation with 
a quantum Bohm-like oc h 2 correction term, for the electrons. The ion continuity and 
momentum equations for the ion fluid density rij and velocity V{ are 

drii d 

— + -(n tUl ) = 0, (3) 
dui dui e dcf) 

l)t +Ui l)x~ = ~ltfdx~ 1 () 
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respectively, where M is the ion mass. Due to their larger mass, the ions are taken as 
non-relativistic and non-quantum. Also, as a first approximation, no thermal effects are 
included for either electrons and ions. The electrostatic potential is obtained from Poisson's 
equation 



d 2 



e 



dx 2 e 



R 2 (dS 



mc 2 \ dt 



+ Hi 



(5) 



where Eo the vacuum electric permittivity, in terms of the appropriate electron charge density, 
as discussed in more detail in Ref. |17j]. Since the interest here is on electrostatic waves 
only, it is sufficient to write the model in terms of 1 + 1 dimensions, with a scalar potential 

= <P(x,t). 

Assuming the electron fluid streaming as a whole through the ionic fluid, an equilibrium 
solution is 

R— J— , S = — / ymc 2 t + p x , rii = n , Ui = , = 0, (6) 
V 7 

in the reference frame where ions are at rest, with the relativistic gamma factor 7 = [1 + 

p 2 / (m^c 2 )] 1 / 2 together with n$ as a the equilibrium density. Notice the modified equilibrium 

electron fluid density, due to the spatial contraction (7 > 1). 

Linearizing the model around the equilibrium ([HD with plane wave perturbations oc 

exp[i(kx — wt)}, one obtains the dielectric function 

1 ^ < [l-h\uj 2 -c 2 k 2 )l{Am 2 c±)} 

u 2 7 [ 7 2 (u;-H 2 -n 2 (u; 2 -c 2 P) 2 /(4m 2 c 4 )] ' { ) 

the dispersion relation being e = 0. In Eq. ([7]), u e = [n e 2 /(me )] 1 / 2 and = 
[rioe 2 1 '{Meq)] 1 / 2 are the electron and ion plasma frequencies, and v = p/{^m) is the beam 
velocity. 

It is instructive to consider the consequences of the dielectric function (J7J) for the three 
separate cases (a) relativistic non-quantum; (b) non-relativistic quantum and (c) joint rela- 
tivistic and quantum. 

III. RELATIVISTIC NON-QUANTUM CASE 

Setting formally h — in Eq. ([7]) yield 

l-F(u,k), (8) 



to 2 u 2 



u 2 7 3 (w — kv ) 2 
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which also defines a non-quantum characteristic function F = F(uj,k). The dispersion 
relation F = 1 is a fourth degree polynomial in the wave frequency u. From the graphics 
of the characteristic function, it follows that stability (four real solutions) is obtained when 
the minimum value -F m i n < 1. When F min > 1 one has two real and two complex conjugate 
roots, one of them an unstable mode. Finally, -F m i n = 1 is marginally stable. In Fig. ([T]) an 
unstable case is depicted. 



















kv 



CO 



FIG. 1: Characteristic function F = F(cj,k) from Eq. (|8|) in the non-quantum case, showing a 
generic unstable equilibrium when F m i n > 1 



Supposing 7(m/M) 1 / 3 <C 1, which in hydrogen plasma is reasonable except in extreme 
relativistic situations (7 > 12.3 or (3 = v/c > 0.997), it is straightforward to obtain 

m \ 1/3 



as the wave frequency for the minimum of the characteristic function. Correspondingly, 

,2 



(9) 



1+7 



M 



> 1 



(10) 



for instability. Hence, relativistic effects (7 > 1) tend to shrink the range of unstable wave 
numbers. 

To obtain the interpretation of the stability analysis, we first rewrite the dispersion rela- 
tion as 

2 ujf(ou-kv) 2 
(uj - UJ n ) u 2 = ^ (11) 



id — UJ r 
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where 

U, P = kV + ^ft, ^ = kV ~^ (12) 

are relativistic versions of the usual Doppler-shifted beam modes. Notice that when ions are 
taken as immobile (M — > oo) one has u PtU as exact normal modes. 

It is interesting to focus on resonant wave numbers and low frequency modes such that 

kv ~ , (13) 

implying uj n ~ 0. One then immediately obtains the estimates 

cu — cj p ~ — 2kv , (u — kv) 2 ~ k 2 v 2 , cu — cu n ~ cu , (14) 

so that Eq. (TiT]) becomes 

" = ~2M^- (15) 
The reason why the dispersion relation (a fourth degree polynomial equation) can be con- 
verted into Eq. fTTS]) . which is of the third degree, is the assumption f[T3"l) . Far from resonant 
or high frequency normal modes cannot be detected in this manner. 
Equation ffTo^) has one real 

\2MJ 77 V ; 

and two complex conjugate roots, 

The plus sign in Eq. (JTTJ) correspond to an exponentially growing mode. Note that the 
growth rate becomes smaller due to relativistic effects. In addition, for all three modes one 
get 



lu ( m 



kv HAff • < 18) 
apart from numerical factors of order unity, so that the low frequency assumption is ac- 
ceptable with the exception of extremely relativistic beams (in the case of hydrogen plasma 
where m/M = 1/1867). Moreover, the instability condition ffTU]) is also fulfilled by the 
resonant wave numbers in Eq. f)13p . 



2l| to be 



The time-averaged energy density <W> in a dielectric medium can be shown 
given by 

<W>= M r (u;)L r , (19) 
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where SE is the electric field perturbation and e r ,oj r are the real parts of the dielectric 
function and wave frequency. In the present dissipation- free problem, e = e r . In Eq. ffT§l) , 
both field and particle energy contributions are already taken into account. 

From the above well-known result, it is found that the energy content of the different 
modes depends on the quantity 



calculated from Eq. ([S]) considering oj = oj r . By inspection, the ion contribution has always 
a positive energy. Concerning electrons, if the wave velocity is larger than the beam's mode 
the energy is positive. From Eq. ( 1121) one has that oj = oj v is a positive energy mode, 
while oj = oj n is a negative energy mode (without loss of generality, we use in this work the 
convention u > 0, k > 0). One can understand the result as follows. When the Doppler- 
shifted frequency of the beam mode is negative, its energy is also negative. Moving through 
the plasma, the beam is then slowed down, losing energy which is the seed for the growing 
amplitude of the wave. 

In conclusion, in this Section the low-frequency instability due to a negative energy 
relativistic beam mode in an oscillating ionic background was investigated. Relativistic 
effects tend to produce a smaller range of unstable wave numbers, as well as a smaller 
growth rate, in comparison to the 7 ~ 1 case. 

IV. NON-RELATIVISTIC QUANTUM CASE 

Setting formally 1/c = in Eq. (J7J) yields 



which also defines a non-relativistic characteristic function F = F(oj,k). The dispersion 
relation F = 1 is again a fourth degree polynomial in the wave frequency oj. However, 
there are qualitative changes in comparison with the non-quantum apparent in Fig. 

(T2J) drawn for v > hk/(2m). It can be verified that the maximum of the characteristic 
function for frequencies kv — hk 2 /(2m) < oj < kv + hk 2 /(2m) is always negative. Hence 
one need F min > 1 for instability, where the corresponding wave frequency oj min satisfy 
< w min < kv — hk 2 / (2m). In passing, the case v < hk/ (2m) can be shown to produce only 
linearly stable oscillations and will be not considered. 



ds 




(20) 




1-F(u,k), 



(21) 
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FIG. 2: Characteristic function F = F(uj, k) from Eq. (|2ip in the non-relativistic case, showing a 
generic stable equilibrium when F m i n < 1. In the graphic, ui± = kv db hk 2 /(2m) 



The dispersion relation F = 1 can be rewritten as 

[u - U) n ) bJ 

where 



, cj [(a; - M 2 - n 2 fc 4 /(4m 2 )] 
cj — w„ 



u„ = kv + \uj 2 p + 



n 2 fc 4 \ 1/2 , / 2 fr 2 A: 4 y /2 
1 , uj n = kv - uj„ H 1 



(22) 



(23) 



are quantum versions of the usual Doppler-shifted beam modes. Notice that when ions are 
taken as immobile (M — > oo) one has u Ptn as exact normal modes. 

In analogy with the previous Section, it is interesting to focus on wave numbers such that 
u n ~ and on low frequencies, or 



2lA\ V2 



Proceeding as before, the result is 



ijj ' 



where 



H 



2M 



0±{H) 



mtr 



measures the strength of quantum effects and where it was defined 



H 



e ± {H) = -j=[i±{i-H 2 Y" 



-1/2 



(24) 



(25) 



(26) 



(27) 



One should compare the relativistic and quantum expressions for the low frequency beam 
modes in Eqs. (1151) and (12"5"j) . respectively. In addition, H is bigger for dense plasmas since 

1/2 

uj e oc n ' . 

The resonance condition can be explicitly solved yielding 

k ^^fy < 28 > 

which gives real wave numbers provided H 2 < 1. In particular, the previous resonant 
wave number is now split in two, due to the quantum recoil term h 2 k A /(4m 2 ) in Eq. (1211) . 
Moreover, one has 



9 + (H) = f (l + f^) + 0(H 5 ) , 0-{E) = 1 - ^ + O(H') , 




(29) 



so that 9± have, respectively, a purely quantum and semiclassic nature. In Fig. ([3D the 
behavior 9±(H) is shown. The two modes coalesce when H = 1. 



FIG. 3: Functions 9± from Eq. (|27[) . Bottom, line: 9+. Upper, dashed: 9- 



Similarly to the non-quantum case, Eq. (125]) can be solved yielding purely oscillatory, 

1/3 

(30) 



fm6±(HY 
w = - — rrV^ I oj p 



damped, 



and exponentially growing 



2M 

1/3 



id 



\ 16M J v ; 



Co' 



'm9±{H) 
16M 



1/3 



Co', 



;i + i\/3) 



(31) 



(32) 



waves. Now one has not only one, but two unstable low frequency modes, according to the 
plus or minus sign chosen in 9±. However, since 9±{H) < 1, it is apparent from Eq. f[3"2j) 
that quantum effects are stabilizing. 

For consistency, it remains to check the low frequency assumption. From Eqs. (l3"0"j) - f[3"iZj) 
it follows that 



kv 



m 



M 



1/3 



^ ~ ( -11 ) eT « i . 



± 



(33) 



Besides, after some algebra one finds from Eq. (I2"T1) . 



m 
M 



9 ± {H) 



1/3 



J? . ~ 
A mm — 



m 



M9 2 ± {H) 



1/3 



> 1 



(34) 



in accordance with the general instability condition. 

Concerning the total energy, one should analyze the quantity 



de 2cu 2 



2u 2 u(u — kv) 



12 • 



(35) 



By using Eq. 
modes, respectively. 



du uj 2 ' {(u-kv) 2 -h 2 ky(4m 2 )Y 
one concludes that oo = u p and uj = u n are positive and negative energy 



V. RELATIVISTIC QUANTUM CASE 

Having in mind the fruitful results following from the low frequency assumption, in this 
work the full relativistic and quantum dispersion relation ([7]) will not be investigated in 
detail. Rather, it is interesting to restrict to slow waves such that 



uj 2 < c 2 k 2 , (36) 



which implies 



. <4 [l + h 2 k 2 /{Am 2 c 2 )] 

uj 2 7 [ 7 2 (w - kv) 2 - h 2 k 4 /(Am 2 )} ' 1 ' 

The qualitative form of the characteristic function derived from Eq. (13TI) is similar to the 
non-relativistic quantum characteristic function from Section IV, with quantitative changes 
due to 1/c 7^ 0,7 > 1. 

Notice that fast wave propagation not satisfying Eq. (1361) can be relevant in some instan- 
ces, but is outside the scope of the present work. For instance, the usual Bohm correction 



is not valid for waves with phase speed greater than the speed of light 22). Moreover, 
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the simplified dielectric function (157)1 does not produce the high frequency pair modes [23] 
admitted by the fully relativistic dispersion relation. 

Once again, the dispersion relation F = 1 is a fourth degree polynomial in the wave 
frequency oj, which can be written as 

w? [(oj - kv) 2 - ^ 2 A;7(4 7 2 m 2 )] 



(oj - 0J n ) u z 



UJ — UJ r 



(38) 



where 



0J r , 



kv + 
kv — 



H 


(■♦ 




+ 


ft 2 fc 4 




4m 2 c 2 y 


47 2 m 2 




(■♦ 


ft 2 *; 2 > 


+ 


n 2 fc 4 




4m 2 c 2 y 


47 2 m 2 



1/2 



1/2 



(39) 
(40) 



Neglecting the ion correction, these are exact beam modes for the dielectric function (1371) . 
They correspond to Doppler-shifted relativistic modifications of the usual Bohm-Pines dis- 
persion relation 24j. 

Similarly to the non-quantum or non-relativistic cases, it is useful to focus on wave 
numbers such that oj n ~ and on low frequencies, or 

1/2 



kv 



1 + 



h 2 k 2 

4m 2 c 2 



h 2 k A 



+ 



47 2 m 2 



oj <^ikv . 



(41) 



The resonance condition gives 



OJ, 



kv 

7 f±(Hn) 



where 



H 

V2 



f3 2 H 2 
47 3 



f3 2 H 2 " 
47 s 



± 



1/2' 



-1/2 



(42) 



(43) 



with (3 = v/c and H ^ given by Eq. 
Proceeding as before, the result is 

,3 



/3 2 7 2 ff 2 
4<4(#, 7 ). 



UJ" 



2 7 4 M 



<p±(h,i) i 



(44) 



For 7 — > l,/3 — > 0, one has Lp±(H ) r y) — > 9±(H). Moreover, the oc /3 2 term in Eq. 
vanishes, so that the quantum non-relativistic results in Eqs. (I25|) and (1281) are recovered. 
On the other hand, choosing (p_(H, 7) (the semi-classical branch) and then letting H — > 
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it is easy to verify that the relativistic, non-quantum results from Section III are recovered, 
considering (p-(H, 7) = 7 5//2 + 0(H 2 ). 

It will be assumed that <p±{H, r y) is not complex, which is true either for 

or for 

27 

P 2 



27 

a 2 



l + /? 2 + ^ + 0(/3 6 ) 



7, 



4 13/3 2 23/3 4 ^ /nB , 



(45) 



(46) 



Equation (j45p extends the condition H < 1 from the preceding Section to the relativistic 
domain, while Eq. (|46p seemingly provides a new parameter regime not existent in the 
non-relativistic approximation. However, in practice it can be verified that the inequality 
(T4"6"|) can be attained only for non-realistic values of H (at least H > 10.8). Hence Eq. (14"5]) 
gives the only important constraint, shown in Fig (J3|) below. We note that the relativistic 
effects allows bigger H values. 




0.5' 



0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 4: The filled area shows the allowable values of the quantum parameter H in terms of (3 = v/c, 
as found from Eq. (|45p 



From Eq. (jHJ) it is obvious that the instability analysis from the preceding Section 
applies, once the replacement 9±(H) — > (±(H,j) is made, where 



C±(#,7) 



r 



(3 2 y 2 H 2 



1 + 



Moreover, generalizing Eq. f[2T?l) one has 



H 



1 + 



(2-7 2 )g 

87 s 



+ 0(tf 5 ), C-(^,7) 



J 



3/2 



1 _ ( 2 ~7 2 )^ 2 

87 s 



(47) 



+ 0(iJ 4 ), 
(48) 
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so that (±(H,j) have, respectively, a purely quantum and semiclassic nature. In Fig. (J5J) 
the behavior of £±(i?>7) is shown for different values of the relativistic parameter 7. The 
two modes coalesce at the maximal quantum parameter found from Eq. fHBT) . Overall, the 
relativistic effects produce a smaller value of (±(H, 7) and hence a smaller instability growth 
rate, which is proportional to [(±(H, 7)] 1 / 3 . 





FIG. 5: Functions C±(H, 7) from Eq. (T3?|) . In all plots, the dashed curve is C- (H, 7), while C+(H, 7) 
is the other one. Parameters are: upper, left, 7 = 1.5; upper, right, 7 = 2.0; bottom, left, 7 = 2.5; 
bottom, right, 7 = 3.0 

Finally, we note that u = u p [Eq. (1311 ] can be shown to be a positive energy mode, while 
uj = u n [Eq. fHU]) ] is a negative energy mode. 

VI. CONCLUSION 

Here we give some notes on the validity of our model. One question in the approach refers 
to the magnetic field generation due to the unbalanced current at equilibrium. Following 

r~l 

the discussion of Bludman et al. [13] we can examine the neglect of static magnetic fields 
according to: 

V x B = —fi n ev B ~ fi n eLv , (49) 
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where B is the magnetic field, fj,Q is the vacuum permeability and L is some characteristic 
dimension of a beam of finite cross section. From this one can estimate the Larmor frequency 
ul as 

eB n e 2 fi Lv u 2 e Lv u L (5Lu e ,_ n . 
lo l = ~ = — =>- — = . (50) 

7m 7m 7cr u e 7c 

For low frequency, the last quantity may be significant, so that the neglect of a return current 
may be unjustified. However, one can always have Wi/w e <C 1 for sufficiently small L. 

Another question is under which ideal experimental conditions the quantum relativistic 
Buneman instability can be investigated, in particular in view of the quantum effects mea- 
sured by the parameter H in Eq. ( I26p . Using S.I. units one has H = 7.25 x KT 20 ^ 72 / /3 2 , 
which attain non-negligible values only for very high densities. For instance, for f3 = 1/10 
(a weakly relativistic beam) one has H > 1/10 only for no > 1.90 x 10 32 m~ 3 . This corres- 



ponds to the parameter regimes of laser-compressed plasma interaction experiments [251.126]. 
Correspondingly, the predicted stabilization effects can be relevant for extremely short wave- 
lengths. 

As a by-product, the present treatment derived the Doppler-shifted relativistic Bohm- 
Pines dispersion relations in Eqs. (1391) and (1401) . The relativistic effects inhibits dispersion, 
with a lower group velocity for a given wavelength, as can be verified. This prediction should 



27|. 



be tested e.g. with the development of multi-Peta-Watt lasers 

In summary, in this work it was developed an unified methodology for the treatment of 
low-frequency beam-driven instabilities in relativistic quantum plasmas. A detailed com- 
parison was made considering different parameter regimes (formally, h = and 1/c 7^ 0; 
h 7^ and 1/c = 0; both h 7^ and 1/c 7^ 0). The corresponding Buneman instabilities 
were shown to be associated to certain beam-plasma negative energy modes. Both quantum 
and relativistic effects cause a stabilizing effect, with the remark that the quantum effects 
produce a splitting of the original Buneman mode into two modes, due to the quantum re- 
coil term as discussed in Section IV. Finally, observe that the model equations have certain 
limitations in particular because of the semiclassical nature (no quantized fields), see 17] 
for more details. Nevertheless, the dispersion relation (j7|) can be equally found using more 
complete, field-theoretic methods [23j. In this context, the positive-negative energy waves 
found here have a more broad applicability, as long as the low frequency assumption is valid. 
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